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THE  INFLUENCE  OF  FLUCTUATIONS  IN  ACOUSTIC  TRAVEL  TIME 
ON  LOCALIZING  THE  POSITION  OF  AN  UNDERWATER  SOURCE 


by 

Tor  Knudsen  and  Angelo  Lombardi 


ABSTRACT 


This  paper  investigates  the  accuracy  of  an  underwater  acoustic  positioning 
system  in  locating  a  sound  source  that  transmits  at  unknown  times  in 
conditions  in  which  the  acoustic  travel  times  between  source  and  receivers 
fluctuate  widely.  The  system  uses  three  omnidirectional  receivers.  The 
accuracy  of  the  systemv depends  very  much  on  the  relative  positions  of  these 
receivers.  With  well+chosen  relative  receiver  dispositions,  the  position 
of  a  sound  source  at  20  km  range  can  be  located  within  ±  150  m,  even  with 
fluctuations  in  travel  time  of  as  high  as  *  100  ms.  < 


INTRODUCTION 


The  need  to  locate  the  position  of  a  sound  source  in  shallow  water  by  means 
of  three  omnidirectional  receivers  led  to  an  experimental  investigation  of 
the  fluctuations  in  travel  time  between  transmission  and  reception  of  the 
sound.  his  showed  that,  in  severe  multipath  conditions,  the  travel  time 
could  vary  by  about  ±  65  ms  between  transmissions.  The  fluctuations  were 
so  rapid  that  two  transmissions  150  s  apart  could  show  a  difference  in 
travel  time  of  about  130  ms  [1].  If  one  knows  when  the  transmission  took 
place,  this  gives  an  error  of  only  about  200  m  in  locating  the  position  of 
the  sound  source.  However,  if  it  is  not  known  when  the  transmission  took 
place,  the  computation  of  the  source's  position  must  be  based  on  the 
differences  in  the  arrival  times  of  the  same  signal  at  three  points  spaced 
some  distance  apart.  The  existence  of  fairly  large  fluctuations  in  the 
travel  times  themselves  means  that  such  fluctuations  will  have  a  large 
influence  on  the  accuracy  of  locating  the  position  of  the  source. 

This  paper  Investigates  the  effect  of  different  magnitudes  of  travel-time 
fluctuation  on  the  accuracy  of  the  computed  position  of  the  sound  source. 
It  also  shows  how  the  effect  of  the  fluctuations  can  be  reduced  by  placing 
the  receivers  in  optimum  relative  positions. 
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1  MATHEMATICAL  MODEL 


In  the  models  the  receivers  are  positioned  as  shown  on  Fig.  1  at  (xj,yi), 
(X2.Y2K  ^x3»y3^‘  The  source  transmits  a  signal  at  time  t  ,  which  Is 
received  at  the  three  receivers  at  times  tj,  t2  and  t3.  We  therefore 
find: 


[ 


(1/c)  /  (x-xi)2  +  (y-yi)2  =  ti-t 
(l/c)  /  (X-X2)2 *  (y-y2)2  ■  t2-t 
(l/c)  / (X-X3)2  +  (y-y3)2  =  t3-t 


(Eq.  1) 


where  c  ,  the  sound  speed,  is  assumed  to  be  constant  and  Independent  of 
depth,  and  the  sound  Is  assumed  to  travel  In  straight  paths.  In  fact  the 
sound  rays  are  bent  in  response  to  the  sound-speed  profile;  this  can  be 
compensated  for  by  modifying  the  sound  speed  used  in  these  calculations. 
From  Eq.  1  we  eliminate  the  unknown  t  and  get 


t2-ti  *  (l/c)  [  /  (X-X2)2  +  (y-y2)2  -  /  (x-xj)2  +  (y-yi )2 ] 
t3-t|  *  (l/c)  [  /  (x-x3)2  +  (y-y3)2  -  /  (x-xi)2  +  (y-yi)2] 


(Eq.  2) 


or,  as  matrix  equation 

Y  =  m(x,y)  ,  (Eq.  3) 

where 


and  m(x,y)  Is  the  right-hand  side  of  Eq.  2. 


X 
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2  OPTIMUM  RELATIVE  RECEIVER  POSITION  [2] 

We  want  to  estimate  the  position  of  the  sound  source 
T 

xk  =  [x,y]t=tk 


from  a  series  of  observations  of  the  differences  in  time  between  the 
acoustic  transmission  and  the  reception  at  the  three  receivers.  The 
mathematical  model  for  this  is  given  in  App.  7A  of  Jazwinsky  [3]: 

Yi  =  M(i )  *(1,k)Xk  +  Vi  ,  i=l . k  ,  (Eq.  4) 

where  {vf}  is  a  white,  but  not  necessarily  gaussian,  vector  sequence  with 

E{  v-f }  =  0  and  E{  v^  vj  }  =  R,  6  ^ j  ,  Ri  >  0 

and  «{1,k)  is  the  state  transition  matrix. 

The  parameter  Xk  may  be  viewed  as  the  state  at  time  tk  of  the  dynamical 
system 


X1+1  =  *(1+1,  1)  Xi 
Yi  =  M(i)  Xi  +  Vi  , 


(Eq.  5) 
(Eq.  6) 


where  Yi  is  the  measurement  at  time  ti 


We  want  to  find  the  estimate  Xk  that  minimizes  the  error  variance 

A  A  ^  a 

E{(Xk-Xk)T  (Xk-Xk)}  =  tr  E{(Xk-Xk)(Xk-Xk)T}  .  (Eq.  7) 


It  can  be  shown  [2]  that  the  linear,  unbiased  minimum  variance  estimate  Xk 
of  Xk  has  a  covariance  matrix  given  by 


E|(Xk-Xkmk-Xk) 


11  •  t: 


(Eq.  8) 


M 


where 


4 


*T/  j  i,  \  uT/ 1  \  o"l 


um  ah  i  \ 
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is  the  information  matrix. 


So  minimizing  the  error  variance  is  equivalent  to  minimizing  t r  i'1  or 

k,l 

maximizing  the  determinant 


In  our  case  the  sound  source  does  not  move,  so  that 

*(i+l,i)  =  *  =  j^1  0  j  =  I 


reduces  to 


l  MT(i )  RT1  M(i )  =  l  w . 
i-1  1  1=1  1 


(Eg.  10) 


so  that  we  shall  maximize 


Wi  |  «  |  Rj  ||  M(i ) 


(Eg.  11) 


Comparing  Eqs.  3  and  6  we  see  that  the  actual  observation  describes  a 
nonlinear  relation  between  the  source’s  position  and  the  measured  time 
difference.  Equation  6  can  be  found  from  Eq.  3  by  linearization. 


3mn ( x *y )  3mi2^x»y) 


ax 


ay 


M  = 


9n>2i(x,y)  8*^22 ( x »y ) 

3x  3y 


Mu  m12 
M21  M22 


(Eq.  12) 


where 


1 


x-x2 


M11  =  7 


c  »  x-X2)2+(y-y2)2  /7x  -xi  )2  +(y-yi )“ 


x-xi _ ) 

l)2+(v-Vl)2  / 


M12  «  t 


y-y2 


y-yi 


c  \n  x-x2)2+(y-y2)2  /  ( x-xi )2  +(y-yi )^ 


4 
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M21  =  i 


X-X3 


x-xi 


c  i  r  ( X-X3 )2  +  (.y-y3 )2  /  (x-xj)2+(y-yi )' 


M22  =  - 


Y-Y3 


y-yi 


C  '  /  (x-x3)2+(y-y3)2  /  { x-xi )2 +(y-yi )2 


Rut  referring  to  Fig,  1  we  now  find  that  M  is  also  given  by 


M  =  I 

c 


cos  02  -  cos  6i 
cos  03  -  cos  ©i 


sin  0  2  -  sin  0j 
sin  03  -  sin  0 1  J 


(Eq.  13) 


where  ©i  ,  02  and  03  are  the  angles  defined  in  Fig.  1. 


If  we  measure  the  arrival  times  with  a  variance  o^.  we  find  that 


R  = 


2a 


2a 


and  thus 


,  1 


3a! 


2  -1 

-1  2 


and  Eq.  11  then  becomes 


Wi  = 


,44 

3a  c 


{ si n (0 1~© 3 )  -  si n (0 2-0 3 )  +  sin(©2-0i)}2  .  (Eq.  14) 


To  maximize  Eq.  14  we  apply  the  Lagrangian  multiplier  method  with  the 
conditional  function 


<t>  (A )  =  Ai-A2+A3  =  0 


(Eq.  15) 


A1 

A2 

a3 


(01  -  63) 
( © 2  "  83) 

(02  -  ©i) 
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and  the  additional  function 
Z  (a  )  =|  Wj  |  +  A<j>  (a  ) 


(Eq.  16) 


Setting  its  partial  derivatives  to  zero  gives  the  following  condition 
A 1  =  A2 

A2  *  a3 
A1  3  A3  • 

Substituting  these  in  Eq.  15  we  find  that 


A 1  =  a2 

gi  ves 

A3  =  0 

a  2  =  a3 

qi  ves 

A 1  =  0 

A 1  =  A3 

gi  ves 

A2  *  2Ai  , 

where  the  last  solution  obviously  is  the  maximum.  Substituting 

A 1  =  A3  =  a2^2  In  Eq.  14  and  maximizing  with  respect  to  A 2  gives 

A  2  =  240° 
and 

Ai  =  a3  =  120°  . 


The  optimum  relative  receiver  positions  are  therefore  in  the  shape  of  a 
triangle  enclosing  the  source  and  such  that  the  angle  subtended  by  each  leg 
of  the  triangle  is  120°,  as  shown  in  Fig.  2. 


y 


FIG.  2  OPTIMUM  RECEIVER  POSITIONS 
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3  SIMULATION 

A  simulation  program  was  developed  to  determine  how  different  receiver 
geometries  would  influence  the  accuracy  in  locating  the  position  of  a  sound 
source.  The  program  computes  the  theoretical  travel  time  between  the 
source  and  the  receivers  and  perturbates  these  delays  with  random  values 
drawn  from  uniform  distributions.  The  position  resulting  from  these 
perturbated  travel  times  is  then  computed,  and  the  process  repeated  with 
another  set  of  random  numbers  in  order  to  plot  the  error  distribution  in 
two  dimensions. 

Assuming  that  the  real  position  of  the  source  is  (xo.yo)  »  Eq.  1  is  used 
to  compute  the  theoretical  travel  times.  The  arrival  times  tj,  t?  and 
*3  are  then  perturbated  according  to 

dl  *  tj  +  s  UC-1.I3 
d2  =  t2  +  s  UC-1,1] 

d3  =  t3  +  s  UC-1,1]  ,  (Eg.  17) 


where  s  is  a  constant  used  for  scaling  and  U[-l,l]  means  a  number  drawn 
from  a  uniform  distribution  of  random  numbers  between  -1  and  1. 

Equation  2  describes  the  relationship  between  the  position  of  the  source 
and  the  measured  travel-time  differences.  Since  it  is  nonlinear,  we  start 
by  linearizing  around  the  actual  position  X0  ,  which  gives  travel -time 
differences  of 


Y-Y0  =  M{ X-X0)  , 

so  that 

X-X0  =  M-^Y-Yq)  =  M-l  d2-<*l-(t2-ti) 

l_d3-di-(t3-ti)  . 

where  M  is  given  by  Eq.  12. 

The  estimated  position  due  to  the  perturbated  travel  time  is 
X  =  M-1(Y-Y0)+X0  . 


(Eq.  18) 
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We  then  compute  the  new  travel  times  using  X  in  Eq.  1  and  relinearize  M 
around  the  new  position  estimate  X  .  Using  Eq.  18  we  then  compute  a  better 
estimate. 

The  iteration  is  ended  when  the  new  and  the  previously  iterated  positions 
are  within  a  chosen  limit,  usually  set  to  5  m.  The  simulation  is  then 
repeated  for  another  set  of  random  perturbations  by  returning  to  Eq.  17. 


4  RANDOM  PERTURBATIONS  OF  THE  TRAVEL  TIMES 

The  computer  program  was  run  on  several  examples  in  order  to  investigate 
how  the  relative  positions  of  the  receivers  influenced  the  error  in 
estimating  the  sound  source's  position.  The  distance  between  the  outermost 
receivers  was  20  km  and  the  sound  source  was  moved  to  different  positions 
within  a  square  of  30  km  *  30  km  dimensions.  Figure  3  shows  the  results 
for  the  three  mean  propagation  delays  of  ±  25  ms,  ±  50  ms  and  ±  100  ms  when 
the  receivers  were  placed  along  the  y  axis  at  10  km  Intervals;  because 
the  errors  are  symmetric  around  a  line  that  is  parallel  to  the  x  axis  and 
crosses  the  y  axis  at  10  km,  only  half  the  field  is  shown. 

It  is  seen  from  the  figure  that  the  uncertainties  in  the  position  of  each 
sound  source  are  ellipses  centred  around  its  real  position,  which  is  marked 
with  a  cross.  The  error  in  the  estimated  position  increases  very  much  when 
the  independent  fluctuations  in  the  time  of  arrival  of  the  signal  at  the 
three  receivers  increase  from  ±  25  ms  to  ±  100  ms.  At  20  km  range  (same 
distance  as  the  receiver  baseline)  the  maximum  error  has  Increased  from 
about  ±  1  km  to  about  ±3.5  km. 

The  situation  is  greatly  improved  when  the  receivers  are  positioned  in  a 
right  angle,  as  shown  in  Fig.  4  for  the  same  three  mean  propagation  delays. 
The  distance  between  the  outermost  receivers  is  again  20  km  and  it  is  seen 
that  even  with  uniform  random  fluctuations  of  ±  100  ms  in  the  time  of 
arrival  the  maximum  error  is  within  ±  500  m  for  a  sound  source  20  km  from 
the  furthest  receiver. 

The  optimum  relative  receiver  disposition  was  found  earlier  to  be  as 
depicted  in  Fig.  2.  The  results  of  placing  the  receivers  in  this 
configuration  are  shown  in  Fig.  5  for  the  same  three  mean  propagation 
delays.  It  is  seen  that  the  errors  within  the  trlanqle  formed  by  the 
receivers  are  very  small  (within  ±  150  m),  even  with  ±  100  ms  time 
fluctuation.  However,  when  the  sound  source  is  placed  outside  the  triangle 
the  errors  start  to  Increase  considerably. 


5  RANDOM  PERTURBATION  OF  SOUND  SPEED 

Normally  one  may  not  know  the  sound  speed  accurately,  or  It  may  vary  along 
the  path  from  the  sound  source  to  the  receiver  and  also  between  the  paths 
to  the  different  receivers.  This  effect  can  be  Interpreted  as  a  range- 
dependent  fluctuation  in  the  acoustic  travel  times  between  source  and 
receivers.  To  give  a  feeling  for  the  problem  we  assume  that  there  is  an 
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uncertainty  in  the  actual  sound  speed,  given  by  (1500  ±  5)  m/s.  This 
results  in  an  uncertainty  of  ±  50  ms  in  travel  time  for  distances  between 
source  and  receivers  that  are  qreater  than  x  ,  qiven  bv 

x  -  Q»Q5  x  1500  x  1505  22 

5 


Since  the  uncertainty  in  the  sound  speed  can  be  interpreted  as  an 
uncertainty  in  acoustic  travel  times,  the  positioning  error  depends,  as 
before,  on  the  relative  position  of  the  receivers. 


CONCLUSION 


This  paper  shows  that  even  with  fluctuation  in  acoustic  travel  times  of  the 
order  of  ±  100  ms  it  is  possible  to  locate  the  position  of  a  sound  source 
at  20  km  range  within  ±  150  m  by  proper  relative  positioning  of  three 
omnidirectional  receivers. 
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RECEIVERS  PLACED  ON  A  LINE. 

Uncertainty  in  position  due  to  fluctuation  in  propagation 
delay  of 

a)  *  25  ms 

b)  *  50  ms 

c)  *  100  ms 
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FIG.  4  RECEIVERS  PLACED  AT  AN  ANGLE. 

Uncertainty  in  position  due  to  fluctuation  in  propagation 
delay  of 

a)  ±  25  ms 

b)  *  50  ms 

c )  *  100  ms 
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FIG.  5  RECEIVERS  PLACED  IN  THE  APICES  OF  A  TRIANGLE. 

Uncertainty  in  position  due  to  fluctuation  in  propagation 
delay  of 

a)  ±  25  ms 

b)  *  50  ms 

c)  ±  100  ms 
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